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The band structure of CasPbO, which possesses a three-dimensional massive Dirac electron at the Fermi 
energy, is investigated in detail. Analysis of the orbital weight distributions on the bands obtained in the 
first-principles calculation reveals that the bands crossing the Fermi energy originate from the three Pb-p 
orbitals and three Ca,-d x 2_ y 2 orbitals. Taking these Pb-p and 0>&r& x i_ y i orbitals as basis wave functions, a 
tight-binding model is constructed. With the appropriate choice of the hopping integrals and the strength 
of the spin-orbit coupling, the constructed model sucessfully captures important features of the band 
structure around the Fermi energy obtained in the first-principles calculation. By applying the suitable 
basis transformation and expanding the matrix elements in the series of the momentum measured from a 
Dirac point, the low-energy effective Hamiltonian of this model is explicitely derived and proved to be a 
Dirac Hamiltonain. The origin of the mass term is also discussed. It is shown that the spin-orbit coupling 
and the orbitals other than Pb-p and Ca,-d x 2_ y 2 orbitals play important roles in making the mass term 
finite. Finally, the surface band structures of CasPbO for several types of surfaces are investigated using 
the constructed tight-binding model. We find that there appear nontrivial surface states that cannot be 
explained as the bulk bands projected on the surface Brillouin zone. The relation to the topological insulator 
is also discussed. 

KEYWORDS: Dirac electron, inverse perovskite, CasPbO, the first-principles calculation, tight-binding model 



a 

O 

o 



> 

o 

o\ 

<N 

o 



1. Introduction 

A low-energy effective Hamiltonian of a material some- 
times becomes a relativistic Dirac Hamiltonian despite 
the fact that electrons in materials are basically de- 
scribed by a nonrelativistic Hamiltonian and the per- 
turbative treatment of the relativistic effects is generally 
legitimated. Such an emergent "Dirac electron" is known 
to host many intriguing phenomena. The well-known and 
well studied realizations of the emergent Dirac electrons 
in materials are graphene 1 3 (two-dimensional, massless) 
and bismuth 4 (three-dimensional, massive). It is also 
known that an organic compound, a-(BEDT-TTF) 2 l3, 
has Dirac electrons in the band structure near the Fermi 
energy. 5 ' 6 Rather new example of emergent Dirac elec- 
trons are surface states of three-dimensional topological 
insulators, 7 which attract great interests and are exten- 
sively studied in these days. Another new example is 
SrMnBi2, which is claimed to have a two-dimensional 
strongly-anisotropic massive Dirac electron. 8 10 A Weyl 
fermion, which has a linear dispersion but is described 
by 2x2 (not 4x4) matrix is also claimed to exist in 
Y2lr2 07, 11 where the time reversal symmetry is broken 
by a magnetic order. 

Recently, we proposed an inverse-perovskite material 
CasPbO and its family as candidates for new materials 
having Dirac electrons. The first-principles calculation 
shows that a three-dimensional massive Dirac electron 
appears in the low-energy band structure of CasPbO. 12 
A Dirac point, which is defined as the center of linear 
dispersion of the Dirac electron, is located on the T-X 
line. Due to the cubic symmetry of this material, exis- 
tence of a Dirac point on the T-X line implies that there 
are six equivalent Dirac points in the Brillouin zone. The 
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most remarkable feature of the Dirac electron in CasPbO 
is that it appears exactly at the Fermi energy, and the 
bands forming the Dirac electron are the only bands that 
cross the Fermi energy. Namely, the Dirac type linear dis- 
persion is the only feature comes up in the vicinity of the 
Fermi energy. 

This point gives a strong contrast to bismuth, which 
is a well-known and well studied example having three- 
dimensional massive Dirac electron in its band structure. 
In bismuth, not only bands with Dirac type linear dis- 
persion, but also a band with parabolic dispersion in- 
tersects the Fermi energy and gives a usual hole Fermi 
surface. 13 ' 14 As a result, one has to dope some carri- 
ers in order to investigate the properties related to the 
Dirac electrons. 15 With this respect, it is important to 
study CasPbO, which has a simpler band structure than 
bismuth, in order to obtain a deeper understanding of 
three-dimensional Dirac electrons in materials. Further- 
more, high symmetry of the crystal structure of CasPbO 
enables us to have an intuitive view on the origin of Dirac 
electron in CasPbO by analyzing the symmetry of wave 
functions for the states relevant to the Dirac electron. 12 
The main topic of this paper is to show the details of such 
an analysis that was briefly reported before. 12 The origin 
of Dirac electron in bismuth is complicated 16 ' 17 and its 
simple and intuitive picture is still missing. Therefore, 
the analysis of CasPbO demonstrated in the following 
will also give hints for the case of bismuth. 

In this paper, a detailed analysis of the band struc- 
ture of CasPbO is presented. First, we construct a tight- 
binding model that captures important features of the 
band structure near the Fermi energy obtained in the 
first-principles calculation. It is assigned that the three 
Pb-p orbitals and three Ca,-d x 2_ y 2 orbitals [Fig. 3(b)] 
should be included in the tight-binding model by an- 
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alyzing the orbital weight distributions on the bands. 
With the proper choice of the hopping parameters and 
the strength of the spin-orbit coupling, the band struc- 
ture obtained in the first-principles calculation is fairly 
well reproduced by the constructed simple tight-binding 
model. The analysis is farther proceeded and we prove 
that the low-energy effective Hamiltonian is really a 
Dirac Hamiltonian by applying a proper basis transfor- 
mation and by expanding the matrix elements with re- 
spect to the momentum measured from one of the Dirac 
points. It is worth noting that not only the Hamiltonian, 
but also the basis wave functions of the low-energy ef- 
fective model are explicitly obtained. These basis wave 
functions play a role of pseudospins of the Dirac model. 
We will also discuss the mass term of the Dirac Hamil- 
tonian in this material, which has been only briefly men- 
tioned in our previous paper. 12 Especially, the roles of 
the spin-orbit coupling and the orbitals other than Pb-p 
and Cdi-d x 2_ y 2 orbitals are explained. The relation be- 
tween CaaPbO and a topological insulator is also dis- 
cussed. It is shown that CasPbO is not a topological in- 
sulator. Finally, the surface band structures of CasPbO 
are studied using the constructed tight-binding model. It 
is found that there exist nontrivial surface bands that are 
nondegenerate and are not explained as the bulk states 
projected on the surface Brillouin zone. 

This paper is organized as follows. Section 2 is used 
to describe details of the method for calculation. In §3, 
the orbital weight distributions on the obtained bands 
are analyzed and a tight-binding model is constructed. 
Section 4 is devoted for analyzing the obtained tight- 
binding model and proving that its low-energy part is 
actually described by a Dirac Hamiltonian. At the end 
of §4, the origin of the mass term is also discussed. In §5, 
we make arguments on the topological properties and the 
surface band structures. The paper is summarized in §6. 

2. First-Principles Calculation 

The band structure calculation is carried out us- 
ing WIEN2k package, 18 in which the full-potential 
augmented-plane-wave method is implemented. The 
crystal parameters required in the calculation are taken 
from the experimental results. 19 Figure 1 shows the 
crystal structure of CasPbO. Perdew-Burke-Ernzerhof 
generalized-gradient approximation 20 is used for the 
exchange-correlation functional. The parameter choice 
of (i^i^max, G max ) = (9.0, 14.0) is used, while radii of 
spherical atomic region for Ca, Pb, and O are chosen to 
be 2.28, 2.50, and 2.28 (a.u.), respectively. 220 momen- 
tum points in the irreducible Brillouin zone, which are 
reduced from 20x20x20 momentum points in the full 
Brillouin zone, are employed in the self-consistent cycle 
in the present calculation. A spin-orbit coupling is taken 
into account within the spherical atomic region of each 
atom via second variational step. 21 

We have verified that changes in the parameters stated 
above (RK max , G max , radii of spherical atomic region, 
and the number of momentum points) do not modify 
the results in this paper, e.g., appearance of the Dirac 
type dispersion in the vicinity of Fermi energy. Although 
the experimental crystal parameters are used in the cal- 




Fig. 1. (Color online) Crystal structure of CasPbO. 



culation, it should be noted that the previous theoretical 
works showed that the optimized lattice constants for 
CasPbO and its family obtained in the first-principles 
calculation agree well with the experimental data, 22,23 
indicating the consistency between theory and experi- 
ments. 

A special care should be paid in the treatment of the 
spin-orbit coupling. Namely, we should be careful in ap- 
plying the second variational step to heavy elements such 
as Pb in which the spin-orbit coupling is expected to be 
strong. 21 However, we think that this is not a serious 
problem in our calculation. One of the reasons is that 
only the state with the total angular momentum J = 3/2 
is important for the Dirac electron in this material as 
will be explained later, while the second variational step 
mainly causes problems for J = 1/2 state. 21 Another rea- 
son is that the band structure does not change (except for 
unimportant points) even if we use the pseudopotential 
method 24 with fully-relativistic pseudopotential for Pb 
atom in which the spin-orbit coupling is treated better 
than in the second- variational step. 

3. Construction of the Tight-Binding Model 

3. 1 Orbital Character of Each Band 

The obtained band structure of CasPbO is presented 
in Fig. 2 together with the orbital weight distributions 
for Ca <i3 Z 2_ r 2, d xy , d x 2_ y 2 1 d zx j yz , and Pb p orbitals. 
As is shown in ref. 12, there appears a Dirac electron 
on the T-X line, and the Dirac point is marked by an 
arrow in Fig. 2(a). The appearance of a Dirac point on 
the T-X line implies that there are six Dirac points in 
the whole Brillouin zone due to the cubic symmetry of 
this material. In other words, Dirac points are found at 
(fco,0,0), (-fe ,0,0), (0,fc ,0), (0,-fe ,0), (0,0, fe ), and 
(0, 0, — fco). Although it is difficult to see in the presented 
energy scale, there exists a very small gap at the Fermi 
energy, and the emerged Dirac electron is actually mas- 
sive with a very small mass. The magnitude of the mass 
gap is about 15 meV. 12 

Before moving on to the discussion on the orbital 
weight distributions, we explain the symmetry of the lo- 
cal environment of each atom. The local environment of 
the Ca atom has a tetragonal symmetry with its tetrag- 
onal axis directing along a line connecting the Ca and 
O atoms (see Fig. 1). The tetragonal-axis directions of 
the three Ca sites in the unit cell are different from each 
other. Then, it is convenient to introduce local coordi- 
nates for each Ca atom, whose definitions are illustrated 
in Fig. 3(a). Using these local coordinates, the Ca-3<i or- 
bitals can be classified into four groups, i.e., d 3z 2_ r 2, d xyi 
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d x 2_ y 2, and d zx / yz , reflecting the local tetragonal sym- 
metry. Note that x : y, and z appearing in the subscripts 
refer to the local coordinates. On the other hand, the 
local environment of the Pb atom has a cubic symmetry 
and all Pb-6p orbitals (p x , p y , and p z ) are classified into 
a single group named as Pb-p. 

The orbital weight distributions of C^-d^ z 2_ r 2 , d xyi 
d x 2_ y 2, and d zx / yz orbitals and Pb-p orbitals are plot- 
ted in Figs. 2(a)-2(e), where the weights are indicated 
by the width of each band. We can see that in the pre- 
sented energy range, the bands below the Fermi energy 
mainly originate from Pb-p orbitals, while the highly 
entangled bands above the Fermi energy mainly origi- 
nate from Ca-3d orbitals. However, note that the top 
of p-bands (bands originate from Pb-6p orbitals) locates 
above the bottom of d-bands (bands originate from Ca- 
3d orbitals) [Figs. 2(c) and 2(e)]. As will be discussed 
later, this overlap between the p- and d-bands is essen- 
tial for the emergence of Dirac electron. Among the Ca- 
3d orbitals, the "center of mass" of Ca-d 32 2_ r 2 orbital 
weights lies at relatively high energy while that of Ca- 
d x 2_ y 2 lies at relatively low energy. This tendency can be 
understood from the crystal field splitting. Specifically, 
the orbitals whose robes are directed to O or Pb (like 
Ca-d 32 2_ r 2) tend to have higher energy than the other 
orbitals since O and Pb are negatively charged in this 
material. On the other hand, the orbitals whose robes 
are directed to direction other than O or Pb (like Ca- 
d x 2_ y 2~) tend to have low energy. Note that C&-d x 2_ y 2 
and Cdi-d zx / yz orbitals have large dispersion and their 
band widths exceed the energy difference of the center of 
masses of each orbital. 

The most important information obtained from 
Figs. 2(a)-2(e) is that two bands forming the Dirac elec- 
tron near the Fermi energy come from C&-d x 2_ y 2 and 
Pb-p orbitals [see Figs. 2(c) and 2(e)]. For Q&-d x 2_ y 2 or- 
bital, there are two reasons for this orbital to come down 
to the Fermi energy: one is its large dispersion and the 
other is that its center of mass lies at relatively low en- 
ergy. On the basis of this observation, we neglect other 
Ca-d orbitals (d 32 2_ r 2, d xyi and d zx j yz ) in constructing 
a model describing the Dirac electron in this material. 




Fig. 3. (Color online) (a) The local coordinates on the three Ca 
sites in a unit cell, (b) Schematic pictures of the basis wave func- 
tions used in the reduced tight-binding model. (c,d) Definitions of 
hopping integrals between the two ^-orbitals (c), between the two 
d-orbitals and between the p- and d-orbitals (d). 



3.2 Tight-Binding Model 

Now, we construct a tight-binding model that de- 
scribes the low-energy band structure of CasPbO. Based 
on the observations in the previous subsection, six or- 
bitals, namely, three Pb-p orbitals (p x , p y: and p z ) and 
three Ca,-d x 2_ y 2 orbitals in the unit cell are taken as ba- 
sis orbitals. The three Card x 2_ y 2 orbitals locate on Cal, 
Ca2, and Ca3 site, respectively. See Figs. 1 and 3(b) 
for the definitions of Cal, Ca2, and Ca3 sites and the 
schematic pictures for the basis orbitals. Note that x, 
y, and z appearing in the subscripts refer to the local 
coordinates as in Figs. 3(a) and 3(b). Furthermore, due 
to the strong spin-orbit coupling on the Pb atom, the 
spin degrees of freedom should be explicitly treated. As 
a result, we use 12 (= 6 x 2) orbitals 

IPxt), \Pyt), \Pzt), \Pxl), \Pyl), \Pzl), (la) 

|dit), \d 2 1), kfet), l<*a>, \d2 |d 3 1)- (lb) 
Using these local orbitals as the basis set, we construct 
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a tight-binding Hamiltonian that is written as 
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where the indices a and a 7 represent p x ,y,z or <ii,2,3, re- 
spectively. Hopping processes described in Figs. 3(c) and 
3(d) are included in the first term of eq. (2) (details are 
explained soon later). The second term represents the 
spin-orbit coupling, which is assumed here to act only on 
the Pb-p orbitals for simplicity. The spin-orbit coupling 
on the Ca-d orbitals is neglected. The Fourier transfor- 
mation of eq. (2) gives 

H = Yl Saa ' ( k ) C ka C ka> , (3) 
k olol' 

where the indices a and a' run through 1 to 12, cor- 
responding to the 12 basis orbitals in the order writ- 
ten in eqs. (la) and (lb). As has been carried out 
in ref. 12, the matrix elements are transformed by 
attaching momentum-dependent phase factors to the 
basis orbitals as \p x , y ,z&) — >* ^ kx+ky+kz ^ 2 \Px,y, z o~), 
\dia) ->> e ikx / 2 \dia), \d 2 a) ->> e ik y/ 2 \d2<r), and \d 3 a) ->> 
e lkz / 2 \ds(j). This is performed for the matrix elements 
to have simple forms. The transformed basis and matrix 
elements are used hereafter. 

Now, the Hamiltonian can be written in a compact 
form as 
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cdd 
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Here, £^ q and A& are 6x6 matrices, while £^ q is a 3x3 

matrix. £% q (q,q f = p or d) represents usual hopping 
processes. 6 in eq. (5) is due to the fact that the non- 
relativistic Hamiltonian has no terms mixing the spin 
components. represents the spin-orbit coupling that 
is caused by relativistic effects, and has matrix elements 
mixing the spin components. In the following subsec- 
tions, we explain details of each matrix in order. 

3.2.1 Spin-Orbit Coupling 

First, A;, contains the contributions from the spin- 
orbit coupling on the Pb-p orbitals. A^ ^ represents the 
on-site L • S coupling for the p-orbitals and is explicitly 
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Next, A^ 1 ' represents the intersite counterpart of the L S 
coupling and is explicitly written as 



(8) 

where c x 

with c a = cos k a (a = x, y, z). Here the intersite coupling 
means the coupling between the nearest neighbor pair 
of two Pb-p orbitals. (Ca-d orbitals are not involved in 
this coupling.) A^ was neglected in ref. 12 since only A<°) 
plays an important role in reproducing the main features 
of the band structure obtained in the first-principles cal- 
culation. However, a close observation of the band struc- 
ture reveals that the energy splittings between the states 
with J = 1/2 and J = 3/2 take different values at the 
T- and at the R-point, which can be explained only if 
we take account of A^. In the following formulation, we 
keep A^ in order to see its effect, although the actual 
magnitude of Ai in eq. (8) is small. 

3.2.2 p-p matrix elements 

Next, £j^ p contains the on-site energy of the Pb-p or- 
bital and hoppings between the two Pb-p orbitals. Con- 
sidering all the hopping processes between the nearest- 
and the next-nearest-neighbor pairs of Pb-p orbitals as 
shown in Fig. 3(c), £% p becomes 



^xx 


e xy 


£xz) 




£yx 


e yy 


e yz 


1, (9) 


^zx 


e zy 


tzz) 





cVV 
C k 



where 

£xx = £ P + 2t v \c x + 2t P 2C x + 4:tpsc y c z + 4tp5C x c x , (10a) 

£yy = ZtplCy + 2t p2 Cy + 4t p3 C Z C X + 4:tp 5 CyCy, (10b) 

^zz £p 

+ 2t pl c z + 2t p2 c z + 4tp 3 c x c y + 4t p5 c z c z , (10c) 

and e X y = —^tp4S x Sy, €y Z = —^tp4SyS z , and e zx = 
—4tp4S z s x with s a = smk a (a = x,y,z). Note that t P 4 
and t p5 were not included in ref. 12 for simplicity. These 
terms improve the agreement between the band structure 
of the tight-binding model and that of the first-principles 
calculation, although the results do not change qualita- 
tively. 



3.2.3 d-d matrix elements 

Next, £% d contains the on-site energy of the Ca-d or- 
bital and hoppings between the two Ca-d orbitals. Con- 
sidering up to the third-nearest-neighbor pairs of Ca-d 
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orbitals as shown in Fig. 3(d), £j* d becomes 



where 
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was not 



included in ref. 12 again, and is included here for com- 
pleteness and to improve the tight-binding model. 




Fig. 4. (Color online) (a) Band structure of the tight-binding 
model (solid lines) and that obtained in the first-principles cal- 
culation (dotted lines). (b,c) Dispersion relations of the tight- 
binding model around a Dirac point. k x and k y are in the region of 
0.06 < k x < 0.22 and -0.08 < k y < 0.08, respectively, while k z is 
(b), and 0.008 (c). Momenta are represented in the unit of 2ir /a 
where a is the cubic lattice constant. 



3.2.4 d-p matrix elements 

Finally, £ pd represents the hybridization between the 
nearest-neighbor pair of the Pb-p and the Ca-d orbitals, 
which is written as 






£x2 


£;r3 







£t/3 




£^2 






where e x2 = -4it dj ,S|C*, e x3 = 4it dp s*cy_, e y3 = 
-4it dp siic?L, e y i = 4it dp S|C|, e zl = -4it dp S|C|, and 
£^2 = 4it dp sjC| with = sin^ (a = x,y,z). 

3.3 Band Structure of the Tight-Binding Model 

Parameters in £^ q and can be determined from 
the maximally localized Wannier function (MLWF) con- 
structing method. 25 27 However, we find that the param- 
eter set obtained using MLWF, which includes rather 
long-ranged hopping processes, does not reproduce the 
band structure near the Fermi energy if the hopping 
integrals are truncated as in Figs. 3(c) and 3(b). This 
is because the long-ranged hopping processes in MLWF 
are important for Pb-p orbitals since their wave func- 
tions have large spreadings, while the hopping integrals 
in Figs. 3(c) and 3(b) are relatively short-ranged. In or- 
der to overcome this problem, we modify some param- 
eters obtained in MLWF to reproduce the band struc- 
ture. This simplification gives no significant problems in 
the following arguments of the appearance of the Dirac 
electron. 

The determined parameters are e p = —1.462, t p i = 
0.290, t p2 = 0.045, t p3 = 0.035, t p4 = 0.0791, t p5 = 0.055, 
e d = 2.146, t d0 = 0.349, t dl = 0.184, t d2 = -0.195, 
t d3 = 0.047, t dp = 0.302, A = 0.330, and Ai = 0.015 (the 
units are eV). Figure 4(a) illustrates the band structure 
in this tight-binding model compared with the results of 
the first-principles calculation. Comparing Fig. 4(a) and 
Figs. 2(c) and 2(e), we can see that the tight-binding 
model well reproduces the bands with large Ca-d x 2_ y 2 
[Fig. 2(c)] and Pb-p [Fig. 2(e)] orbital weights. This is 
natural since we have taken account of these orbitals in 
constructing the tight-binding model. Note that not only 



the global band structure, but also the most important 
feature of this system, i.e., the Dirac type linear disper- 
sion in the vicinity of the Fermi energy, is captured. 

In order to check this point, Figs. 4(b) and 4(c) pro- 
vide three-dimensional plots of the dispersion relations of 
the tight-binding model around one of the Dirac points, 
(&o,0,0). Figure 4(b) for k z = shows that two cone- 
shaped bands touch with each other at an isolated point 
in the k x -k y plane, while Fig. 4(c) shows that a gap 
emerges between the two bands as k z becomes finite. 
These observations confirm the existence of a three- 
dimensional Dirac electron in this model. However, we 
find that the Dirac electron in this tight-binding model 
is massless while the first-principles calculation gives 
a small mass gap of 15 meV. We find that this mass 
term originates from the contribution from the other or- 
bitals which are not included in the present tight-binding 
model. This will be discussed later in §4.4. 

4. Low-Energy Effective Hamiltonian 

In this section, we concentrate on the Dirac point on 
the £^-axis, namely, the Dirac point at (0,0, /co). This is 
because we are familiar with choosing the z-axis as the 
quantization axis of the spin. When we discuss the band 
structure around other Dirac points, such as (&o,0,0), it 
is convenient to direct the spin quantization axis along 
the x-axis. 

4.1 Analysis of the Tight-Binding Model 

In order to prove that the low-energy effective Hamil- 
tonian of this material is really a Dirac Hamiltonian, 
we first introduce a new basis set. Here, the eigenstates 
at the T-point are taken as the new basis set since the 
Dirac point under consideration is relatively close to the 
T-point. In addition, the high symmetry of the T-point 
allows us to classify the eigenstates easily. Explicitly, the 
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Fig. 5. (Color online) Weight distributions of the new basis 
states on the bands along the fc z -axis. (a)-(f) show the weights 
of k>3 3), \p3_1), \pi 1), \d\ a), \dLa), and \di,a), respectively. 
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new basis wave functions are defined as 
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Here, \pjj z ) is just the wave function having a total angu- 
lar momentum J with its z-component J z realized from 
the on-site spin-orbit coupling of the Pb-p orbitals. The 
bases \d[a), \d 2 a), and \d 3 a) originate from the Ca-d or- 
bitals [eq. (15)], and spin up and spin down components 
are not mixed since the spin-orbit coupling on the Ca 
atoms is neglected in this model. 

Figure 5 shows the weight distribution of the new basis 
states on the bands along the k z -ax.is. From this figure, it 
is clearly seen that the bands crossing the Fermi energy 
and forming the Dirac electron come from the \p33_) and 
\d' 2 <j) states. Therefore, it will be enough to keep only the 
four states |p3_3_), |p3i), \d' 2 t) 5 and \d 2 \) out of 12 states 



in eqs. (14) and (15) in deriving the low-energy effective 
Hamiltonian. When we examine Fig. 5 more closely, we 
can see that a small weight of \d' 3 a) state is mixed in 
the band forming the Dirac electron [Fig. 5(f)]. However, 
the mixed weight is so small that it gives no harm in the 
following arguments. 

4.2 Expansion around Dirac point 

Reflecting the arguments in the previous subsection, 
we derive the low-energy effective Hamiltonian in the 
following two-step process. In the first step, the matrix 
elements of £% q are expanded with respect to k x and 
k y up to the first order (remember that we are concen- 
trating on the Dirac point on the /c^-axis), and at the 
same time, the basis set in eq. (1) is transformed into 
that in eqs. (14) and (15). In the second step, matrix 
elements irrelevant to |P| 3), \p33_), \d 2 t), and \d 2 \) are 
dropped. This procedure is slightly different from that 
used in ref. 12, but the two procedures are essentially 
the same. In the following, we make use of the notation 

k-\- == kg; zL iky . 

In the first step, i.e., expanding the matrix elements 
with respect to k x and k y (or fc±), and changing the basis 
set, £k is transformed as 



c'PP 



gfdd 



Here, t'™ 



can be explicitly written as 



(16) 



c'PP 



9i 


95k- 








95k- 





9bk+ 


92 








94 


-g'lk- 








92 


-95k- 




-94 








-95K 


91 







95 k + 


94 







93 








-95K 


-94 


g',k- 





9s 



with 



51 = A + 2A! + / , 

52 = A + |Ai(l + 2c z ) + i(/ 



2/i), 



53 = -2A - |Ai(2 + c z ) + i(2/ + /i), 



7 

(17) 

(18a) 
(18b) 

(18c) 



g 5 = 4t P 4S z /a/3, 
and g' b =g 6 /y/2, g'l = \/3g' 6 , and 
fo = (e p + 2t pl ) + (2t p2 + 4t p5 )(l + c z ) 
h = (cp + tt p2 + 4^) + (2t pl + 8t p5 )c 
On the other hand, £' k dd can be explicitly written as 
(g% 0\ 



(18d) 
(18e) 



(19a) 
(19b) 



ctdd 










97 99 

99 9s 

g 6 













97 99 

99 9s/ 



(20) 
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with 



96 — h ~ h •> 



^ = -(/3 + 2/ 4 + /5-4/ 6 ), 



08 



09 



3(2/3 H 
^(f 

-g-(/3 



/ 4 + 2/ 5 +4/ 6 ), 
fi + h - h)i 



and 



h = + 2t d i + 2£ d2 ) + 2t d2 c 2 , 
/ 4 = (e d + 4t d2 ) + 2t d ic z , 
/s =4(^ + 2t d 3C z ), 
/e = 4(t d0 + 2t d3 )c j • 
The off-diagonal matrix, S'^ p , is expressed as 
(—h\k- hik + h%k + 



(21a) 
(21b) 

(21c) 

(21d) 

(22a) 
(22b) 
(22c) 
(22d) 



\ 



-h[k+ 






-ft' 2 fc_ 







-fc' 2 'fc_ 



K 



/i 2 ^+ h^k 



3^+ 


-h 3 k- 

-Kk + 



J 



-h[k 

-h A 
-h x k+ -h 2 k_ 
-ti A h'lk- -h' 2 'k+ 



(23) 

with hi = itdpCj, ft 2 = i^ p (2 + C|)/V / 3, /13 = — 2it rfp (l — 
c§ )/v / 6, and ft 4 = 8itd P s*/>/3, where ft' lj2j3 = ^1,2,3/^, 

In the second step, we keep only the matrix elements 
related to bn), |d 2 t)> and 1^2 Di and drop all 

the other matrix elements. As a result, the Hamiltonian 
is transformed to a 4 x 4 matrix which can be written 

as 

( 9l h 2 k+ \ 
gi -h 2 k 
91 
\ -ft5fe+ # 7 J 

Finally, suppose that g\ = gj is satisfied at some k 
and expand gi and gj as 



(24) 



e 



(25) 



where = k z — ko and eo = gi (= 07) at k z = fco- 
Substituting eq. (25) into eq. (24), the Hamiltonian be- 
comes, 



-CpSk z 







c p Sk z 



h* 2 k+ 



-h 2 k+ 


c d Sk z 




\ 



(e + £c£fe)i + 



I -cSk z 


V 





—cSk z 




Cd5k z J 

-h 2 k+ \ 

h 2 k. 

cSk z 

C(5fe 2 J 



(26) 



where c and Sc are defined as c = + c p )/2 and c = 
— c p )/2. At this point, it is easy to demonstrate that 
eq. (26) is really a (tilted) massless Dirac Hamiltonian 



by changing the order of the rows and columns. Here, 
we emphasize that the four basis states of this Dirac 
Hamiltonian come from \ps 3), \ps 3), \d' 2 t)> an d |cZ 2 i)- 

In deriving eq. (24), we have neglected the states \p3 1 ) , 
|j9 3 1 ) , |pi 1 ), |_p 1 1), and |rf]_cr) . Here, let us discuss why 
we can neglect these states in detail. Firstly, the energies 
of these states are away from the energy of the Dirac 
point (eo). For \p\ 1) and \p k 1), this is mainly because 
1 03 -e | is large at fe^ = fe^ . For \pz 1), |jP|i), and l^cr), 
this is due to the term /14 in eq. (23). Namely, although 
\g 2 —€Q I and |^6 — e o| a t fc^ = k z o are not so large, the term 
/i4, which is finite even for |fc±| =0, makes eigenenergies 
for |_p| 1 ) , |_p 3 1 ) ? and |c?icr) away from eo- Secondly, we 
find that the matrix elements between these states and 
the states that we have kept in eq. (23) are at most in 
the first order in k±. This means that the contributions 
of these states to the low-energy Hamiltonian are at least 
in the second order in k±. 

On the other hand, we should be careful in neglect- 
ing Idgcr), since the matrix elements between \d 3 cr) and 
\d 2 cr} are finite even for \k±\ =0 [see eq. (23)]. In order 
to eliminate these constant off-diagonal matrix elements, 
we have to choose an appropriate linear combination of 
\d 2 cr) and \d 3 a) instead of pure \d 2 cr) as a basis state. 
However, using such a linear combination only results 
in a renormalization of /i 2 ,3, h' 2 3 , and h 23 in eq. (23) 
since the transformation leading to the linear combina- 
tion of \d 2 a) and \d 3 a) only mixes the second and the 
third columns (or the fifth and the sixth columns) of the 
matrix in eq. (23). As a result, the low-energy effective 
Hamiltonian should have the same form as eq. (24). Fur- 
thermore, the weight of \d 3 a) mixed in the band form- 
ing the Dirac electron is very small as we have seen in 
Fig. 5(f). 

4.3 Origin of the Dirac Electron 

As we discussed in ref. 12, the origin of the Dirac elec- 
tron in this model can be clarified by considering the 
effects of the hybridization between the p- and d-bands 
(tdp) on the band structure. As we have pointed out, 
the top of p-bands locates above the bottom of d-bands. 
Then, if the hybridization td p is neglected, the p- and 
d-bands should form usual hole- and electron-like Fermi 
surfaces around the T-point, and there appear no Dirac 
electrons. Although the p- and d-bands cross, the cross- 
ing points are not isolated points in the Brillouin zone. 
Then, we consider the effect of finite td p . In general, finite 
tdp causes the band anticrossing of the p- and d-bands, 
and eliminates the crossing points between the p- and 
d-bands. However, in the present model, the band anti- 
crossing does not occur on the T-X line because of the 
symmetry of wave functions (ref. 12 and see below). Con- 
sequently, the crossing points are left on the special and 
isolated points in the Brillouin zone that are located on 
the r-X line. This means that there emerge Dirac points 
in this model. Note that these arguments have some sim- 
ilarities with those in ref. 28. 

Finally, we address why the band anticrossing does not 
occur on the T-X line. For simplicity, we concentrate on 
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the Dirac point on the /c^-axis, which is one of the T-X 
line. In this case, the states relevant to the Dirac electron 
are Iph), and \d' 2 <j). As defined in eqs. (14) and 

(15), and \p$ 3) are the superpositions of \p x cr} and 

\Py(r), while \d 2 a) is a superposition of \d\<j), l^c), and 
\cl3a). The latter wave functions schematically depicted 
in Figs. 6(a) and 6(b). From Fig. 6, it is apparent that 
\p x a) (or \p y cr)) cannot hybridize with the wave func- 
tions of Figs. 6(a) and 6(b) by a symmetrical reason. 
This means that there are no matrix elements between 
|_p 3 3) (or |^3 1)) and \d 2 a) on the /c^-axis, which explains 
the absence of the band anticrossing. Note that if a mo- 
mentum is off from the /c^-axis, there appear finite matrix 
elements between (or |p3_i)) and \d 2 cr). 



— t 


V 1 - 


+ 


+ 

'(a) 



377 



J 



(b) 





Px 1 


1 Py 






d Z x | 


J d JZ 






d x 2_ y 2 


(c) 





Fig. 6. (Color online) (a,b) The wave functions involved in \d' 2 a) 
states for the case that the momentum is on the fc^-axis. (c) 
Schematic pictures for the orbitals on the Pb site. 



4.4 Origin of the Mass Term 

The above arguments nicely explain the origin of the 
massless Dirac electron in the tight-binding model. How- 
ever, the first-principles calculation shows that there ex- 
ists a small but finite mass gap at the Dirac point in 
CasPbO. In this subsection, we clarify the origin of this 
mass term. First, let us check the irreducible represen- 
tations to which the bands on the T-X line belong. As 
shown in Fig. 7, it turns out that the two bands forming 
the Dirac electron belong to the same irreducible repre- 
sentation Tj. If they belong to different representations, 
it is impossible to have a mass gap. However, since they 
belong to the same representation, finite matrix elements 
between these two bands can exist, which result in a fi- 
nite mass gap. 




Fig. 7. (Color online) Irreducible representations of the bands 
on the T-X line. The notations used in ref. 29 are used to label 
irreducible representations. 



that the simple tight-binding model introduced in the 
previous subsection is a good starting point discussing 
the mass term. In order to include the small mass term 
to the effective model, some minor effects should be in- 
cluded. In the following, we consider the effects caused 
by the orbitals not included in the previous subsection. 
As a representative example, we add the Pb-bd orbitals 
to the basis set. (The inclusion of other orbitals leads to 
the similar conclusion.) Since these orbitals are located 
far away from the Fermi energy (about —15 eV below 
the Fermi energy according to the first-principles calcu- 
lation), we treat these orbitals in a perturbative way. 
We concentrate on the effective model around the Dirac 
point on the A^-axis as before in the subsequent argu- 
ments. 

The Pb-bd orbitals considered here are also shown in 
Fig. 6(c). From Fig. 6(c), we can see that Pb-d zx (Pb- 
d yz ) can hybridize with \p x a) (\p y cr)) in the z-direction 
if k z is finite. As a result, and |pss), which are 

relevant to the Dirac electron, can hybridize with Pb-d zx 
and Pb-d yz , since they are composed of \p x a) and \p y (j). 
On the other hand, we can see that Pb-d x 2_ y 2 orbital can 
hybridize with the wave functions depicted in Figs. 6(a) 
and 6(b). As a result, \d 2 <j), which is the other state 
relevant to the Dirac electron and is composed of wave 
functions in Figs. 6(a) and 6(b), can hybridize with Pb- 
d x 2_ y 2 orbital. To summarize, and \d' 2 o) 

will be modified in a perturbative way as 



|p§§> — > N§) + aT |pw i > 



|PaI> — > |P3 3> + -r i |Pb-d I ) 



\d' 2 a) — > \d' 2 a) + -^\Pb-d x 2_ y 2<j) , 



Wl 

Ai 

W2 

A" 2 



(27) 
(28) 



(29) 



where 



|Pb-ds> 



t>+i|Pb-d^t>), (30) 



\Pb-d 



i> = ^(|PW. 



i) - i\Ph-d yz i) 



(31) 



Here, wi^ and Ai^ represent the matrix elements and 
the level differences between the originally included and 
newly added orbitals, respectively. Note that w\ and W2 
depend on k z , and wi is zero at k z = while W2 is finite 
at k z = from the symmetry of the orbitals. 

Now, assuming the standard on-site L ■ S coupling for 
Pb 5d orbitals (not Pb 6p orbitals), there appear new 
matrix elements between \ps 3), and \d' 2 o) orbitals. 

Explicitly, we find that 

(4 t I^Lsbff ) ~ a(Pb-d x 2_ y 2 t |^ L s|Pb-d f ) = 0, 

(32a) 
(32b) 

. Similarly, we have (d 2 | |^ls|P||) ~ 
\^2a\ and (d 2 | |^ls|Ps 3) = 0. Then, adding these con- 



(4 t |^ls|P3 3) - a(Pb-d x 2_ y 2 t |^Ls|Pb-^3 

2 2 



where a = 



Since the mass term is very small in CaaPbO, we think 



J. Phys. Soc. Jpn. 



FULL PAPERS 



tributions to eq. (26), the effective Hamiftonian becomes 



4" = (e + ScSk z )l + 



( 
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-c5k z 


-h£k- 
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—cSk z 
m 



-h 2 k + 

m 
cSk z 




m 
h 2 k. 


cSk. 



\ 



J 



( 33 ) 

with m ~ y/2a\. This is nothing but a (tilted) massive 
Dirac Hamiltonian. 

We can obtain the two important features about the 
mass term from the above derivation. First, the mass 
term is small since it requires the inclusion of the orbitals 
locating far away from the Fermi energy. Second, the 
mass term scales with the coefficient for the spin-orbit 
coupling, A. 

We must note that any other orbitals having the same 
symmetry as Pb-d zx , d yzi d x i_ y i orbitals can be sources 
of the mass term. For example, some states originated 
from the Ca-3d orbitals (other than d x 2_ y 2, which is in- 
cluded in the original basis set) have symmetries of Pb- 
d zxi and d yz orbitals. Thus, inclusion of the Ca-3d or- 
bitals other than d x 2_ y 2 will also contribute to the mass 
term. Again, these orbitals are located away from the 
Fermi energy, and thus, the mass term induced by these 
orbitals are also small. The resultant mass term scales 
with A also in this case. 

5. Topological Property and the Surface Band 
Structure 

5.1 Topological Property 

In this subsection, we discuss the relation of CasPbO 
to the topological insulator. As we have stated, the top of 
p-bands, whose wave functions have odd parity at the Y- 
point, locates above the bottom of d-bands, whose wave 
functions have even parity at the T-point. This kind of 
structure, i.e., an overlap of the two bands with opposite 
parity, is called an "inverted" band structure, which can 
leads to a topological insulator in some cases. 30 ' 31 Fur- 
thermore, Pb atom has a very strong spin-orbit coupling. 
Then, it is tempting to assign CasPbO as a topological 
insulator. In fact, CasPbO and its relatives are suggested 
as potential topological insulators in ref. 32. However, 
unfortunately, CasPbO is not a topological insulator as 
explained in the following. 

The Z 2 topological number v$ for CasPbO is readily 
evaluated since this material has an inversion symmetry. 
According to ref. 30, vq can be obtained from the formula 



8 N 



(-ir=n*=n n^(r<), 



(34) 



2=1 m=l 



for the materials having an inversion symmetry. In this 
formula, Ti represents the time reversal invariant mo- 
menta (TRIM) in the Brillouin zone and £2™^) repre- 
sents the parity of the wave function for the 2m-th state 
at r^. TV is the number of filled bands counted neglecting 
the Kramers degeneracy. For the simple cubic crystal of 
Ca 3 PbO, TRIM are the T-, X-, M-, and R-points where 
there are three X- and M-points. Using WIEN2k package, 
we derive the irreducible representations and the parities 
of the states at TRIM. Here, we concentrate on the states 
whose energies are in between —4 eV and eV, since the 



bands below —4 eV are associated with the completely 
filled shells of the involved atoms such as 0-2p orbitals 
and give only trivial contributions to i/q. There are three 
bands between —4 eV and eV, and thus, N = 3. The 
obtained results are summarized in Table I. From the 
table, we can see that the right hand side of eq. (34) is 
equal to 1, which implies that this material is topologi- 
cally trivial. 



Table I. Irreducible representations and the parities of the states 
whose energies are in between — 4 eV and eV at TRIM. 
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We can understand why CasPbO is topologically triv- 
ial as follows. As we have checked, there exists an "in- 
verted" band structure at the T-point. However, the in- 
verted bands at the T-point have a double degeneracy 
other than Kramers degeneracy. In other words, they 
have a quadruple degeneracy if the Kramers degeneracy 
is included. As a result, this "inverted" band structure 
gives ( — l) 2 = 1 in the right hand side of eq. (34). This 
makes CasPbO topologically trivial. In order to turn this 
material into a topological insulator, at least the degen- 
eracy at the T-point should be lifted by lowering the sym- 
metry of the crystal. 

5.2 Surface Band Structure 

Although we have proved that CasPbO is not a topo- 
logical insulator, we study the surface band structure of 
CasPbO. The surface band structures are calculated us- 
ing the tight-binding model in a slablike lattice structure 
having two (top and bottom) surfaces. Here, surfaces are 
treated as simple terminations of the lattice points for 
simplicity (i.e., we use open boundary conditions). In the 
following, we consider the cases in which the slab has 001 
or 111 surface. In CasPbO, there are two types of 001 sur- 
faces. The schematic pictures for these two 001 surfaces 
are shown in Fig. 8(c) and 8(d). These two types of sur- 
faces are named as 001 A [Fig. 8(c)] and 001B [Fig. 8(d)], 
respectively. As we can see from Figs. 8(c) and 8(d), 001 A 
surface has equal number of Ca and Pb sites, while 001B 
surface has only Ca sites. (In Fig. 8(d), O sites are also 
written, but O sites are not included in the basis set of 
the tight-binding model.) In contrast, there is only one 
type of 111 surface [Fig. 9(b)]. Actually, there is another 
possible 111 surface, but that surface has only O sites 
and is equivalent to the 111 surface in Fig. 9(b). 

The surface band structures for the 001 surfaces are 
summarized in Figs. 8(a)-8(c). The top and bottom sur- 
faces of the slab used in the calculation are 001 A for 
Fig. 8(a), 001B for Fig. 8(b), and 001 A and 001B for 
Fig. 8(c). Figure 8(f) shows the surface Brillouin zone 
for the 001 surface and the arrows indicate the path on 
which the band structures in Figs. 8(a)-8(c) are shown. 
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Fig. 8. (Color online) Surface band structures for the 001 sur- 
faces. The top and bottom surfaces of the slab are (a) 001A, (b) 
001B, and (c) 001A and 001B. (d,e) Schematic pictures for (d) 
001A and (e) 001B surfaces, (f) Surface Brillouin zone for the 001 
surfaces. Dots represent the Dirac points projected to the surface 
Brillouin zone. Allows represent the path on which the band struc- 
ture is plotted. 



The positions of the Dirac points projected on to the sur- 
face Brillouin zone are represented by dots. Considering 
the relation between the arrows and dots in Fig. 8(f), 
the shapes of the bulk continuum (the region filled by a 
large number of bands) in Figs. 8(a)-8(c) are easily un- 
derstood: the bulk continuum has a gap except near the 
projected Dirac points. 

Very interestingly, we can see that several bands ap- 
pear near (0,0) in all of Figs. 8(a)-8(c) apart from the 
bulk continuum. We can see that these separated bands 
are definitely surface states since they are affected by 
the choice of the surfaces. When we closely investigate 
the surface states, we find that they are doubly degener- 
ate in Figs. 8(a) and 8(b). This degeneracy is a Kramers 
degeneracy since the slab used in Figs. 8(a) and 8(b) 
has an inversion center because the top and bottom sur- 
faces are identical. On the other hand, the surface states 
in Fig. 8(c) are nondegenerate. This is because the slab 
does not have an inversion center since the top and bot- 
tom surfaces are different for Fig. 8(c). From these re- 
sults, we can see that the nondegenerate surface states 
exist when we look at one of the top or bottom surface. 
In this sense, this model resembles to a topological insu- 
lator, although we have checked that CasPbO is not a 
topological insulator. Note that it is natural for the sur- 
face bands to have no Kramers degeneracy (if we look 
at one surface), since the spin-orbit coupling is included 



in our model and the existence of a surface breaks the 
inversion symmetry. 

The surface band structure for 111 surface is shown 
in Fig. 9(a). Again, we can see the surface states in 
Fig. 9(a). As in the case of the 001 surface, the surface 
bands are nondegenerate if we look at one of the top or 
bottom surface. These results suggest that the nonde- 
generate surface states in this model are rather robust, 
although this system is not a topological insulator. It is 
a very interesting future problem to clarify the relation 
between the topological insulators and CasPbO series. 




Fig. 9. (Color online) (a) Surface band structure for the 111 sur- 
face, (b) Schematic picture for the 111 surface, (c) Surface Brillouin 
zone for the 111 surface. Dots and allows represent the same things 
as in Fig. 8. 



6. Summary and Discussion 

In this paper, we first constructed a tight-binding 
model that describes the low-energy band structure of 
CasPbO. The three Pb-p orbitals and three G^-d x i_ y i 
orbitals are included in the basis set, reflecting the re- 
sult of the analysis of the orbital weight distributions 
on the bands. The appropriate choices of the hopping 
integrals lead to the satisfactory agreement between the 
band structure obtained in the first-principles calculation 
and that of the constructed tight-binding model. Espe- 
cially, the Dirac electrons, which is the most important 
feature of the band structure of CasPbO, emerge in the 
tight-binding model. 

It was then confirmed that the low-energy effective 
Hamiltonian in the vicinity of the Fermi energy is really a 
Dirac Hamiltonian by applying the basis transformation 
and the expansion of the matrix elements with respect to 
k± . In this procedure, four basis wave functions required 
for the realization of the Dirac Hamiltonian are identified 
as |ps s), \ps 3), \d' 2 t)j an d \d' 2 I). Each basis wave func- 
tion corresponds to each components of the pseudospin of 
the Dirac electron. Therefore the information of the basis 
wave functions clarified in this paper will play important 
roles in considering the possible scheme for controlling 
the internal degrees of freedom of the emergent Dirac 
electrons. 

We have also made a detailed argument on the origin 
of the mass term of the Dirac electron. The conclusion 
is that the orbitals other than Pb-p and Ca-d x 2_ y 2 or- 
bitals and the spin-orbit coupling play key roles in pro- 
viding the mass term. From this analysis, we can easily 
understand the reason for the smallness of the mass term. 
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Namely, the finite mass term requires the inclusion of the 
orbitals whose weights are lying far away from the Fermi 
energy. Since the spin-orbit coupling on the Pb atom is 
essential for the mass term, the mass term can be con- 
trolled by replacing Pb atoms by, for example, Sn. Specif- 
ically, the alloy of Ca3(Pbi_ ;E Sn ;E )0 will give a smaller 
mass gap since the spin-orbit coupling of Sn is weaker 
than that of Pb. In fact, the first-principles calculation 
for CasSnO results in the smaller mass gap compared 
with that obtained in CaaPbO. 12 

Finally, we have discussed the surface band structures 
of CaaPbO using the constructed tight-binding model. 
It is found that there exist nontrivial surface bands that 
cannot be explained as the bulk states projected on the 
surface Brillouin zone. Note that ref. 11 shows that there 
should be nontrivial surface states for materials having 
Weyl fermions. Interestingly, the observed surface states 
are nondegenerate if we look at one surface. Further stud- 
ies on the surface band structure in the present model, 
such as to investigate the relation between CasPbO and 
topological insulators, are highly desired. In this paper 
surfaces were treated as simple open boundaries. When 
we compare the present results and some future experi- 
mental results, it will be required to consider the effect 
of surface potential since the surfaces in Figs. 8(c), 8(d), 
and 9(b) are not charge neutral. 

The successful construction of the simple tight-binding 
model for CasPbO encourages further search for new ma- 
terials having Dirac electrons. The detailed analysis pro- 
vided in this paper are useful in the future experimental 
and theoretical studies of the CasPbO and related ma- 
terials. 
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